Numerical prediction of blast‐induced stress wave from large‐scale underground explosion     

Chengqing Wu  Yong Lu  Hong Hao 《国际地质力学数值与分析法杂志》2004,28(1):93-109

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

非线性Rossby波所满足的三阶Zakharov方程与低频振荡     

罗德海 《成都信息工程学院学报》1993,(4)

本文在基本气流具有水平切变的情况下,利用摄动法导出了非线性Rossby波所满足的三阶Zakharov方程,然后,考虑了基流具有弱切变的情况,通过使用三阶Zakharov方程研究了Rossby波列的第一类不稳定性问题。结果表明:通过非线性作用,大气中的Rossby波列可产生调制不稳定。同时,本文对这种不稳定的区域,增长率和周期进行了详细的计算,并讨论了波振幅、波数、纬度和基流切变对它们的影响,指出Rossby波列的调制不稳定可以激发30～60天的低频振荡。  相似文献   

Loss of water from Mars:: Implications for the oxidation of the soil     

H. Lammer  H.I.M. Lichtenegger  I. Ribas  E.F. Guinan  S.J. Bauer 《Icarus》2003,165(1):9-25

The evolution of the martian atmosphere with regard to its H₂O inventory is influenced by thermal loss processes of H, H₂, nonthermal atmospheric loss processes of H⁺, H₂⁺, O, O⁺, CO₂, and O₂⁺ into space, as well as by chemical weathering of the surface soil. The evolution of thermal and nonthermal escape processes depend on the history of the intensity of the solar XUV radiation and the solar wind density. Thus, we use actual data from the observation of solar proxies with different ages from the Sun in Time program for reconstructing the Sun's radiation and particle environment from the present to 3.5 Gyr ago. The correlation between mass loss and X-ray surface flux of solar proxies follows a power law relationship, which indicates a solar wind density up to 1000 times higher at the beginning of the Sun's main sequence lifetime. For the study of various atmospheric escape processes we used a gas dynamic test particle model for the estimation of the pick up ion loss rates and considered pick up ion sputtering, as well as dissociative recombination. The loss of H₂O from Mars over the last 3.5 Gyr was estimated to be equivalent to a global martian H₂O ocean with a depth of about 12 m, which is smaller than the values reported by previous studies. If ion momentum transport, a process studied in detail by Mars Express is significant on Mars, the water loss may be enhanced by a factor of about 2. In our investigation we found that the sum of thermal and nonthermal atmospheric loss rates of H and all nonthermal escape processes of O to space are not compatible with a ratio of 2:1, and is currently close to about 20:1. Escape to space cannot therefore be the only sink for oxygen on Mars. Our results suggest that the missing oxygen (needed for the validation of the 2:1 ratio between H and O) can be explained by the incorporation into the martian surface by chemical weathering processes since the onset of intense oxidation about 2 Gyr ago. Based on the evolution of the atmosphere-surface-interaction on Mars, an overall global surface sink of about 2×10⁴² oxygen particles in the regolith can be expected. Because of the intense oxidation of inorganic matter, this process may have led to the formation of considerable amounts of sulfates and ferric oxides on Mars. To model this effect we consider several factors: (1) the amount of incorporated oxygen, (2) the inorganic composition of the martian soil and (3) meteoritic gardening. We show that the oxygen incorporation has also implications for the oxidant extinction depth, which is an important parameter to determine required sampling depths on Mars aimed at finding putative organic material. We found that the oxidant extinction depth is expected to lie in a range between 2 and 5 m for global mean values.  相似文献   

Monte-Carlo simulation of Mercury's exosphere     

Peter Wurz  Helmut Lammer 《Icarus》2003,164(1):1-13

Because of its proximity to the Sun and its small size, Mercury has not been able to retain its atmosphere and only a thin exosphere surrounds the planet. The exospheric pressure at the planetary surface is approximately 10⁻¹⁰ mbar, set by the Mariner 10 occultation experiment. The existence of gaseous species H, He, and O has been established by Mariner 10. In addition Na, K, and Ca have been observed by ground based instrumentation. Other elements are expected to be found in Mercury's exosphere since the total pressure of the known species is almost two orders of magnitude less than the exospheric pressure.It is intended to measure these exospheric particle densities in situ with an instrument on board of ESA's BepiColombo Mercury Planetary Orbiter (MPO) spacecraft. Since the expected exospheric densities are very small we developed a Monte-Carlo computer model to investigate if such a measurement is feasible along the MPO spacecraft orbit. We model energy and ejection angle distributions of the particles at the surface, with the emission process determining the actual distribution functions. Our model follows the trajectory of each particle by numerical integration until the particle hits Mercury's surface again or escapes from the calculation domain. Using a large set of these trajectories bulk parameters of the exospheric gas are derived, e.g., particle densities for various atomic and molecular species. Our study suggests that a mass spectrometric measurement is feasible and, at least at MPO's periherm, all species that are released from the surface will be observed.  相似文献   

A close one-term approximation to the highest Stokes wave on deep water     

R.C.T. Rainey  Michael S. Longuet-Higgins   《Ocean Engineering》2006,33(14-15):2012-2024

  相似文献   

Titan's atmospheric engine: an overview     

Robert E. Samuelson 《Planetary and Space Science》2003,51(2):127-145

Heating occurs in Titan's stratosphere from the absorption of incident solar radiation by methane and aerosols. About 10% of the incident sunlight reaches Titan's surface and causes heating there. Thermal radiation redistributes heat within the atmosphere and cools to space. The resulting vertical temperature profile is stable against convection and a state of radiative equilibrium is established. Equating theoretical and observed temperature profiles enables an empirical determination of the vertical distribution of thermal opacity. A uniformly mixed aerosol is responsible for most of the opacity in the stratosphere, whereas collision-induced absorption of gases is the main contributor in the troposphere. Occasional clouds are observed in the troposphere in spite of the large degrees of methane supersaturation found there. Photochemistry converts CH₄ and N₂ into more complex hydrocarbons and nitriles in the stratosphere and above. Thin ice clouds of trace organics are formed in the winter and early spring polar regions of the lower stratosphere. Precipitating ice particles serve as condensation sites for supersaturated methane vapor in the troposphere below, resulting in lowered methane degrees of supersaturation in the polar regions. Latitudinal variations of stratospheric temperature are seasonal, and lag instantaneous response to solar irradiation by about one season for two reasons: (1) an actual instantaneous thermal response to a latitudinal distribution of absorbing gases, themselves out of phase with the sun by about one season, and (2) a sluggish dynamical response of the stratosphere to the latitudinal transport of angular momentum, induced by radiative heating and cooling. Mean vertical abundances of stratospheric organics and aerosols are determined primarily by atmospheric chemistry and condensation, whereas latitudinal distributions are more influenced by meridional circulations. In addition to preferential scavenging by precipitating ice particles from above, the polar depletion of supersaturated methane results from periodic scavenging by short-lived tropospheric clouds, coupled with the steady poleward march of the continuously drying atmosphere due to meridional transport.  相似文献   

Seismic rays in an ellipsoidal earth model: theoretical analysis and estimation of deviations due to the ellipticity     

Jean-Francois Guéguen  Pho Hoang-Trong  Hilaire Legros 《Geophysical Journal International》1987,89(3):839-855

Summary. A first-order form of the Euler's equations for rays in an ellipsoidal model of the Earth is obtained. The conditions affecting the velocity law for a monotonic increase, with respect to the arc length, in the angular distance to the epicentre, and in the angle of incidence, are the same in the ellipsoidal and spherical models. It is therefore possible to trace rays and to compute travel times directly in an ellipsoidal earth as in the spherical model. Thus comparison with the rays of the same coordinates in a spherical earth provides an estimate of the various deviations of these rays due to the Earth's flattening, and the corresponding travel-time differences, for mantle P -waves and for shallow earthquakes. All these deviations are functions both of the latitude and of the epicentral distance. The difference in the distance to the Earth's centre at points with the same geocentric latitude on rays in the ellipsoidal and in the spherical model may reach several kilometres. Directly related to the deformation of the isovelocity surfaces, this difference is the only cause of significant perturbation in travel times. Other differences, such as that corresponding to the ray torsion, are of the first order in ellipticity, and may exceed 1 km. They induce only small differences in travel time, less than 0.01s. Thus, we show that the ellipticity correction obtained by Jeffreys (1935) and Bullen (1937) by a perturbational method can be recovered by a direct evaluation of the travel times in an ellipsoidal model of the Earth. Moreover, as stated by Dziewonski & Gilbert (1976), we verify the non-dependence of this correction on the choice of the velocity law.  相似文献   

A numerical investigation of the destruction of peroxy radical by Cu ion catalysed reactions on atmospheric particles     

Howard B. Ross  Kevin J. Noone 《Journal of Atmospheric Chemistry》1991,12(2):121-136

To determine if Cu mediated reactions on atmospheric particles are important to HO₂ chemistry in the ambient atmosphere, Cu molalities were calculated from measured Cu aerosol particle concentrations, mass and number size distribution data from a site in central Sweden. A comparison of characteristic times indicates that at low relative humidities the reaction is limited by the mass transport of gas phase HO₂ to the particle surface and not by the chemical kinetics of the aqueous reaction. Comparison of half-lives for particle reactions and the gas phase destruction of HO₂ to form H₂O₂ indicate that heterogeneous reactions on aerosol particles may have important consequences on the chemistry of HO₂ and H₂O₂ in the troposphere.  相似文献   

密度和压缩系数的散射层析成像法   总被引：1，自引：1，他引：0  

袁晓晖  王妙月 《地球物理学报》1991,34(6):753-761

本文在速度成像的基础上研究了同时对密度和压缩系数成像的散射波层析成像法.对不同散射角度的计算可以得到一系列反演图像,拟合这些图像,从而可以有效地达到对密度和压缩系数（或速度）成像的目的.与单纯的速度成像相比,增加了反演的难度.首先是对资料的方位性要求增加;其次是对资料的利用率下降.即便如此,从对较少量的炮点和检波点资料的数值计算来看,仍取得了满意的成像结果.我们对组成字母“A”的散射体结构进行了成像计算,结果能够同时再现密度和压缩系数,成像清晰,表明了方法的可行性,并能应用于复杂结构的成像问题.  相似文献   

青藏高原及邻近区域的S波三维速度结构   总被引：25，自引：5，他引：20  

周兵  秦建业 《地球物理学报》1991,34(4):426-441,T001

本文收集了WWSSN台网和我国台网中13个地震台站的长周期地震记录,用140条10-90s瑞利波频散曲线和作者提出的Tarantola-Backus面波频散层析成象方法,作了青藏高原及邻区的速度反演,得出该地区岩石层速度结构的三维图象.结果表明,1.在10-110km深度范围内,速度结构出现与大地构造特征相关的分区性,显示出四个构造单元:青藏块体、柴达木-巴颜喀拉-三江块体、塔里木块体和印度块体.2.高原内部,深度为10-70km内速度较低,莫霍界面呈不对称盆形分布,藏北那曲附近地壳厚度超过70km,高原边缘壳厚为45-50km,90-110km为高速异常,表明高原内部存在上地幔盖层.3.高原北部的班公湖断裂和东部的三江断裂系是该区重要的分界线,是岩石层结构存在明显差异的重要接触部位,可能是冈瓦纳古陆与欧亚古陆的缝合带.4.柴达木-巴颜喀拉-三江块体内部速度分布不均匀,地壳厚度由北向南从45km加深到60km;在深度90-110km存在一低速层.5.塔里木地块内速度随深度均匀增加,从地壳到上地幔110km内没有发现低速层.地壳厚度约50km.  相似文献